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We report a hierarchical first-principles investigation on the entangled effects of lattice dimensionality and 
bond characteristics in the lattice dynamics of silicene and germanene. It is found that bond bending (stretch¬ 
ing) negatively (positively) contributes to Griineisen constant y, which results in the negative acoustic (positive 
optical) y. The layer thickening (bond weakening) caused by chemical functionalization tends to increase (de¬ 
crease) the acoustic (optical) y, due to the increased (decreased) bond-stretching effect. The excitation of the 
negative-y modes results in negative thermal expansion, while mode excitation and thermal expansion compete 
with each other in thermomechanics. The sensitive structural and electronic responses of silicene and germanene 
to functionalization help us to derive a generic physical picture for two-dimensional lattice dynamics. 

PACS numbers: 68.35.bg, 67.63.-r, 63.22.-m, 65.40.De 


Two-dimensional silicon and germanium (silicene and ger¬ 
manene) llh have attracted great interest after the advent of 
graphene 101 . The low dimensionality of graphene |0 ren¬ 
ders its charge density, electronic states, and various physical 
and chemical properties readily controllable | 00 ], which is 
highly helpful for advanced materials and devices. Apart from 
having these low-dimensional benefits, silicene (Si) and ger¬ 
manene (Ge) are also expected to be easily incorporated into 
existing silicon-based industry nil. 

Si and Ge have buckled hexagonal lattices 13, 0, where 
quantum Hall effects, valley polarization, tunable band gap, 
and fast intrinsic mobility have been discovered IMl. This 
indicates their promising applications in electronics, optics, 
and fundamental-physics research. Si has been synthesized 
on various substrates (Ag CMl, Ir Cl, and ZrBz 0), 
while Ge is still under exploration. Monolayer Si transistors 
operating at room temperature have been fabricated recently 
ylfl. Si and Ge can be easily functionalized in various envi¬ 
ronments. Chemical functionalizations (e.g., hydrogenation) 
can controllably open a ba nd g ap in Si and Ge by changing the 
orbital hybridization 12214251] . and hydrogenated germanene 
(HGe) with a band gap of 1.53 eV has been synthesized ll2l . 
The supporting substrates also have similar functionalization 
effects on Si and Ge 11271431|] . These fast synthesis and fab¬ 
rication progresses lead to an urgent demand for the knowl¬ 
edge of the thermal expansion and thermomechanics of pris¬ 
tine/functionalized Si and Ge, because the accumulated ther¬ 
mal strain and stress may influence the performance and life¬ 
time of devices working at finite temperatures. 

In this letter, the phonon spectra, Griineisen constants, ther¬ 
mal expansion, and temperature-dependent stiffness of Si and 
Ge are investigated by first-principles simulation. The roles of 
low dimensionality, layer thickness, and bond strength are dis¬ 
entangled by analyzing the effects of chemical functionaliza¬ 
tion. The revealed physical picture is useful for understanding 
the lattice dynamics of various two-dimensional systems. 

The electronic structure and phonons are calculated us¬ 
ing density-functional methods 13211 . The temperature depen¬ 
dent lattice constant (a) is calculated using the self-consistent 


quasiharmonic approximation (SC-QHA) method 
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where k and A are the wavevector and branch indices for a 
phonon mode; Nk is the k-point number; Uk,A, and jkp 
are the total electronic energy, and the internal energy and 
Griineisen constant (y - of a phonon mode, respec¬ 

tively. To guarantee that y is calculated between two modes 
with the same symmetry, the phononic k-p theory 03411 is used 
to sort the phonon branches. Eq. O is self-consistently 
solved, and the phonon frequencies are updated after each iter¬ 
ation step (Aw = -y^ca). The thermal-expansion coefficient 
(a - f ^) can be re-expressed as 
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where Cy is the isovolume heat capacity and 820 is the two- 
dimensional bulk modulus 
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The contributions of thermal expansion and phonon excita¬ 
tion to B 2 D are included in (ii - ^ 2 ) and (53 + 54 ), respectively, 
and the phonon excitation is omitted in the quasistatic bulk 
modulus { 8 * 20 ) 13511 . Temperatures higher than 600 K are not 
considered here, due to the low thermodynamic/kinetic sta¬ 
bility of these systems found in experiments |i7, 2^ 13], as 
well as to the appearance of high-order anharmonicity that not 
included in QHA 137, 
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Hydrogenated Si and Ge (HSi and HGe) are used to study 
the effects of chemical functionalization. The structures are 
shown in Fig. [11^)^ the calculated equilibrium lattice 
constants (a), buckling heights (Aj), and bond lengths {d) are 
listed in Tab. [B which are consistent with other experimental 
lUllllSIllS and theoretical results ISIHEl . The dif¬ 
ferential electron densities between the pre- and post-bonding 
states (Ap) are projected onto the z direction (Ap^), which are 
shown in Fig. [Hb) and (c). The cr and tt bonds coexist and 
compete with each other in buckled Si and Ge due to the par¬ 
tial sp^ orbital hybridization iBi. Therefore, apart from 
the space between the Si/Ge atoms (|z| < A^), the electrons 
also have some observable accumulation (Ap^ > 0) outside 
(|z| > Aj.) due to the n bonding (Fig. mb, c), left panels). 
Si has a smaller Aj, (0.22 A) than Ge (0.34 A), indicating the 
less sp^ hybridization in the former. Hydrogenation height¬ 
ens Aj (Tab. nil and removes the n electrons (Fig. [Hb, c), right 
panels), due to the increased sp^ hybridization. The bond en¬ 
ergy {eb) is defined to be the energy cost to break a Si-Si/Ge- 
Ge bond, which is used to indicate the bond strength. Al¬ 
though the cr bond is strengthened by hydrogenation, the total 
strength (e*) of a Si-Si/Ge-Ge bond is decreased (Tab. |I]i, due 
to the removal of the tt bond. On the other hand, hydrogena¬ 
tion increases the dynamical stability of the bond, due to the 
removal of the ct-tt competition. These changes in the bond 
characteristics (hybridization, geometry, strength, and dynam¬ 
ical stability) by hydrogenation will have profound effects on 
the lattice dynamics. 


A detailed knowledge on the vibrational states is a prereq¬ 
uisite for understanding various phononic and thermodynamic 
properties. The phonon dispersions and density of states (gpi,) 
of Si, HSi, Ge, and HGe are shown in the left and middle pan¬ 
els of Fig. |2] where the phonon branches are labeled according 
to their initial symmetries at the F point, and the amplitudes 
of the in-plane (XY) and out-of-plane (Z) vibrations are in¬ 
dicated by the line widths. The evolution of the phononic 
eigenvectors of Si along the FK path are visualized in Fig. 
|3] The TA, TO, and LO branches in Si/Ge only consist of 


TABLE I: Equilibrium lattice constant (a), zero-point expansion 
{S.p), buckling bight (A.), bond lengths (4x-x and dx-u, X =Si or 
Ge), and bond energy (eh) of Si, HSi, Ge, and HGe. 



«(A) 

S,p (%) 

A. (A) 

1 

>< 

dx-H (A) 

eh (eV) 

Si 

3.87 

0.13 

0.224 

2.28 

- 

2.91 

HSi 

3.89 

0.11 

0.359 

2.36 

1.501 

2.51 

Ge 

4.04 

0.12 

0.340 

2.43 

- 

2.41 

HGe 

4.07 

0.13 

0.367 

2.46 

1.552 

2.00 


the XY vibration, and the ZA branch of the Z vibration (Fig. 
Ha, c)). Thus, the atomic displacements in these four modes 
(Fig.0 are only modulated by the Bloch phase factor (e 
However, in the LA and ZO branches, the XY and Z vibra¬ 
tions hybridize and swap with each other at k ~ 0.22FK/FM 
in Si (0.40FK/FM in Ge). This vibrational hybridization is 



FIG. 1; (a) Structures and (b, c) Ap^ of pristine/hydrogenated Si and 
Ge. The origin (cr and/or n bondings) of the electron accumulation 
(Apj > 0) near the Si/Ge atoms are indicated in shades. 


caused by the nonorthogonal covalent bonds, which also has 
been observed in functionalized graphene ll3^ and M 0 S 2 |34| 
with nonorthogonal C-C and Mo-S bonds, respectively. In 
planar graphene, the XY modes (TA, LA, TO, and LO) are 
completely decoupled with the Z modes (ZA and ZO) S, 
because the C-C bonds in the XY plane are orthogonal to the 
Z direction. The LA-ZO vibrational hybridization changes 
the LA mode from XY to Z symmetry, which renders the 
LA branch able to couple with the ZA branch at the Brillouin 
zone boundary (Fig. |2a, c)), where they are similar in both 
frequency and symmetry. The LA-ZO and LA-ZA couplings 
result in the flat dispersion of the LA branch at 110 cm ' in 
Si (60 cm ' in Ge). The acoustic/optical modes with Z vi¬ 
bration have lower frequencies than their counterparts with 
XY vibration, because the former ones mainly consist of the 
out-of-plane bending of the Si-Si/Ge-Ge bonds, which are 
softer than the later ones consisting of the in-plane bond bend¬ 
ing/stretching. This is a low-dimensional effect in phonon 
spectra. Ge has a larger buckling height Aj than Si, and the 
contribution of bond stretching in the ZO branch is larger in 
more buckled layer, which tends to stiffen the ZO modes. This 
is a thickness effect in the lattice dynamics of low-dimensional 
systems 03311 . Therefore, the LA-ZO coupling in Ge is 
smaller than that in Si, resulting in less XY (more Z) vibration 
in the LA (ZO) branch and larger dispersive LA branch in Ge 
(Fig. I2c)). Hydrogenation additionally introduces four bend¬ 
ing (TB, TB*, LB, and LB*) and two stretching modes (ZS and 
ZS*), and the H atoms also have some resonating displace¬ 
ments in the acoustic and optical modes (Fig. |3b, d)). Hy¬ 
drogenation increases the layer thickness (A^, Tab. |I|l, which 
stiffens the ZO modes and enlarges the LA-ZO frequency gap. 

ratio is increased from 0.51 (0.51) in Si (Ge) 


The -Jdmz: 

((Jzo+(^la)/2 

up to 0.69 (0.65) in HSi (HGe), which significantly decreases 
the LA-ZO coupling, and the increases XY (Z) vibration in the 
LA (ZO) branch. With less Z vibration in the LA branch af- 
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FIG. 2: Phonon dispersions, gpi,, and 7 dispersions of (a) Si, (b) HSi, (c) Ge, and (d) HGe. The line widths in the phonon dispersions are scaled 
by the amplitudes of XY and Z vibrations, and the antisymmetric H modes in (b) and (d) are labeled with the superscript 



r^. 



FIG. 3: The phonon modes in Si at different k points. 


ter hydrogenation, the LA-ZA coupling at the zone boundary 
is then decreased due to their enlarged symmetry difference, 
making the LA branch there more dispersive. In addition, the 
ZA branches in both Ge and HGe have wiggles near the T 
point, which indicates their low dynamical stability (originat¬ 
ing from weak cr bond and cr-rr competition). This agrees with 
the facts that Ge still has not be synthesized due to the pre¬ 
ferred germanium-substrate alloying ill, and that HGe has 
a low amorphodization/decomposition temperature (<348 K) 
i^ . HGe is synthesized prior to Ge, which may be due to the 
dynamical-stability enhancement by hydrogenation. 

The Griineisen-constant ( 7 ) dispersions are shown in the 
right panels of Fig. |2l 7 (ZA) has large negative values, which 
is ubiquitous for the bending ZA modes in two-dimensional 
systems lf3l . These bending ZA modes resemble the flexu¬ 


ral vibration in guitar string, where a tensile makes the lat¬ 
tice/string stiffer to such bending vibration, resulting in a 
higher frequenc y an d then a negative 7 . This is called guitar¬ 
string analogy MTOl] or membrane effect M- Apart from 
7 (ZA), 7 (TA) and 7 (LA) also have negative values in Si and 
Ge. In a TA mode, the lattice is transversely distorted in the 
XY plane (Fig. O, and the restoring force comes from both the 
in-plane bending and stretching of the covalent bonds. Bond 
bending negatively contributes to 7 due to the guitar-string 
analogy, while bond stretching has a normal positive contri¬ 
bution to 7 . When the bond is not strong enough, the effect 
of bond stretching can not compete with that of bond bend¬ 
ing, resulting in negative net y(TA). Diamond and graphene 
only have positive 7 (TA) 11391 kill due to the strong C-C 


bonds, while bulk silicon and germanium have negative 7 (TA) 
142. 43 I] due to the relatively weak Si-Si and Ge-Ge bonds. 
Through the LA-ZO and LA-ZA couplings, the weight of 
Z vibration in the LA branch increases when leaving the F 
point (Fig. Ilia, c), left panels), which increases the contribu¬ 
tion of out-of-plane bond bending, and then results in nega¬ 
tive 7 (LA), especially at the zone boundary. All the optical 
branches only have positive 7 , due to the dominating effect of 
bond stretching, which is similar in functionalized graphene 
SEIl]. However, the 7 (ZO) in planar graphene is negative, 
because of the dominating bond-bending effect in those ZO 
modes. Except for the rippling part of the 7 (ZA) dispersion in 
Ge, hydrogenation increases the 7 of ZA, TA, and LA modes 
in Si and Ge, because the layer thickening increases the pos¬ 
itive contribution of bond stretching in these acoustic modes. 
The abnormal rippling part of the 7 (ZA) in Ge is caused by 
the low bond dynamical stability, and those ripples in 7 (ZA) 
are largely suppressed by hydrogenation, because the related 
ZA modes become closer to normal two-dimensional bending 
modes. The 7 (ZO) of Si is exceptionally high (= 9) near the 
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r point, due to the stretching of the nonorthogonal strong n 
bonds in the lowly-buckled Si lattice (Tab. |I|i. Hydrogenation 
more or less decreases the j of ZO, TO, and LO modes in 
Si and Ge, because the elimination of the n bond decreases 
the bond strength and then the effect of bond stretching. The 
bending and stretching hydrogen modes (TB, TB*, LB, LB*, 
ZS, and ZS*) in HSi and HGe all have near-zero y, because 
the restoring forces on those perpendicular C-H bonds are in¬ 
sensitive to the lateral-size variation. This means that hydro¬ 
genation is an ideal approach to study the effect of chemical 
functionalization on the lattice dynamics of Si and Ge, be¬ 
cause those modes specific to H atoms have negligible effects 
on the lattice anharmonicity. 

The above analyzed phononic properties can help under¬ 
stand various thermodynamic properties of Si and Ge, as well 
as the effects of functionalization. The calculated thermal- 
expansion coefficient (a), isovolume heat capacity (Cy), and 
quasiharmonic/quasistatic bulk modulus ( 820 !Bhj) of Si, HSi, 
Ge, and HGe are shown in Fig. |4] A fully excited phonon 
branch contributes 1.0 ks to Cy, and the nominal number of 
excited branches is about three at 100 (50) K in both Si and 
HSi (Ge and HGe) (Fig. IH^b)). This means that the excita¬ 
tion of the three acoustic branches (ZA, TA, and LA) dom¬ 
inates at temperatures lower than 100 (50) K, and the opti¬ 
cal modes (ZO, TO, and LO) are only considerably excited at 
higher temperatures. Acoustic modes mostly have negative y, 
while optical modes only have positive y (Fig. |2] right pan¬ 
els). Due to the subsequent excitation of these negative- and 
positive-y modes, a firstly decreases with increasing tempera¬ 
ture until 100 (50) K (Fig. EJa)), and then increases. However, 
a is still consistently negative, due to the dominating contri¬ 
bution from the negative-y modes. Although the excitation 
of the additional H modes makes the Cy of HSi (HGe) larger 
than that of Si (Ge) above 120 K (Fig. EJb)), its effect on the 
minimum-a temperature is negligible due to their near-zero y 
of these modes (Fig. |2b, d), right panels), and hydrogena¬ 
tion also has a negligible effect on the expansion caused by 
the zero-point vibrations (d^p, Tab. |I|i. Therefore, the hydro¬ 
genation effect on the lattice anharmonicity only comes from 
the chemical functionalization on the Si-Si/Ge-Ge bond. Al¬ 
though the hydrogenation-induced decreases of 820 (by 30%, 
Fig. lUc)) and optical y (Fig. |2a, b)) tend to further decrease 
the negative a (Eq. |2]), Q'(HSi) is still close to a(Si), due to the 
canceling effect from the increase of acoustic y (Fig. |2la, b)). 
Qr(Ge) is much higher than Q;(HGe) (by 4 x 10'®K * above 
300 K), because y(ZA) is significantly decreased by hydro¬ 
genation (Fig. |2c, d)), which negatively contributes to a. 

Thermal expansion is caused by the excitation of anhar- 
monic phonons, while both thermal expansion and phononic 
excitation have individual contributions to the thermomechan¬ 
ics (Eq. |3] and |4|i. Both contributions are considered in the 
temperature dependence of 820 , while only the former in that 
of 82 P,. The negative thermal expansion results in the in¬ 
crease of 8212 under heating, due to the increased curvature 
of the potential-energy surface under lattice contraction (i.e., 

dB* 

< 0). The stiffening effect from the zero-point vibrations 



T (K) T (K) 



T (K) T (K) 


FIG. 4: Temperature dependence of (a) a, (b) Cy, and (c, d) 620 and 
ffjo of Si, HSi, Ge, and HGe. In (c) and (d), the curve slopes (in 10“^ 
eVA^^K^*) at 300 K are indicated. 


makes 820 larger than at low temperatures. However, 
B 2 £)(T) curves have lower increasing rates than their 8 * 2 j^(T) 
counterparts, and the 820 of Si even decreases under heat¬ 
ing, which is due to the softening effect from the excitation 
of negative-y modes (54 term in Eq. |3]l. Therefore, thermal 
expansion and phonon excitation have reverse effects on 82 D 7 
which also has been observed in positively expanding M 0 S 2 
03411 . Although hydrogenation decreases 820 and 
the bond weakening (Tab. |I|l, hydrogenation increases both 


the slopes of B 2 d(T) and B 2 q(T) curves and -^). 
equals a-^ x a, and the of Si at 300 K is increased by 
5.0x10“^ eVA^^K '. As Si and HSi have the same a, this in¬ 
crease is only ascribed to the enlargement of potential-surface 
anharmonicity, namely, to the decrease in (from -6.3 to 


-13.3 eVA ^). The increase in the of Si by hydrogenation 

is even larger (by 6.7x10“^ eVA^^K^') than due to the 


increased acoustic y (^4 ~ 2 Cy^yicj). The decrease in the 
negative a of Ge by hydrogenation contributes to the increase 
of ^ (by 3.1x10-^ eVA-^K-' at 300 K), which is partially 
canceled by the increase in (from -18.6 to -13.8 eVA“^). 

The increase in ^4 (by 2.4x10“^ eVA^^K^') is smaller than 

dB* 

that of due to the decrease in y(ZA) after the dynamical- 
stability enhancement by hydrogenation (Eig. |2lc, d)). 


In summary, the roles of lattice dimensionality and bond 
characteristics in the lattice dynamics of silicene and ger- 
manene, as well as in the chemical functionalization effects, 
have been revealed. Eor the substrate-supported silicene and 
germanene 1127143111 , as well as other two-dimensional mate¬ 
rials, the environmental interactions may bring some compli¬ 
cated structural and electronic influences 11441 - 1461] . The physi¬ 
cal picture and analysis methods established here for the cor- 
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relation between structure, electronics, and lattice dynamics 
will be useful for studying the corresponding lattice anhar- 
monicity therein. 
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